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Abstract 

Motivated by the regime of QCD explored nowadays at LHC, where both the total 
energy of collision and momenta transfers are high, we investigate evolution equations of 
high energy factorization. In order to study such effects like parton saturation in final 
states one is inevitably led to investigate how to combine physics of the BK and CCFM 
evolution equations. In this paper we obtain a new exclusive form of the BK equation 
which suggests a possible form of the nonlinear extension of the CCFM equation. 

1 Introduction 

The Large Hadron Collider (LHC) is already operational and Quantum Chromodynamics 
(QCD) is the basic theory which is used to set up the initial conditions for the collisions 
at the LHC as well as to calculate hadronic observables. The application of perturbative QCD 
relies on the so called factorization theorems which allow to decompose a given process into a 
long distance part, called parton density, and a short distance part, called matrix element. Here 
we will focus on high energy factorization [HE]. The evolution equations of high energy fac- 
torization sum up logarithms of energy accompanied by a strong coupling constant, i.e. terms 
proportional to a™ ln m s/sq, which applies when the total energy of a scattering process is much 
bigger than any other hard scale involved in a process. 

Until now, in principle, the BFKL [3H5], BK [6HE] and CCFM PEHE] evolution equations 
were used on equal footing since the energy ranges did not allow to discriminate between these 
frameworks. However, there were indications already at HERA [13] for the need to account for 
nonlinear effects in gluon density. These observation are supported by recent results obtained 
in [Ff|IT5]. On top of this, the results from [16] point at the need to use the framework 



which incorporates hardness of the collision into BFKL like description. With the LHC one 
entered into a region of phase space where both the energy and momentum transfers are high. 
Therefore, one should provide a framework where both dense systems and hard processes can 
be studied. This might be achieved with a relatively simple nonlinear extension of the CCFM 
equation where one could take into account both the gluon production and recombination in 
the description of final states. 

In this paper we obtain such an equation, see eq. (122]) . In order to arrive at this equation 
in the first crucial step we perform resummation of virtual and unresolved real contributions 
in the BK equation in a similar fashion as in [17]. As a result, we arrive at a new form of the 
BK see eq. Qlfip equation where both the linear and nonlinear parts are folded with a Regge 
form factor in which singularities of the linear part have been resumed. In this new form, the 
singularity in the unresolved real contribution of the s-channel real gluon is canceled by the 
virtual contribution. This is the minimal condition in order to eventually perform a Monte 
Carlo simulation based on the BK equation. 

However, the BK equation concerns inclusive observables and does not allow for applica- 
tions to the description of exclusive final states. Thus, we need the CCFM equation which is 
applicable to the description of the exclusive processes. We propose the nonlinear extension 
of the CCFM equation being motivated by the resumed form of the BK equation in which we 
replace the Regge form factor by the non-Sudakov form factor and introduce angular ordering 
(coherence). In addition we supplement the BK kernel with a large x part. For a different 
approach to the extension of the CCFM equation to allow for gluon saturation we refer the 
reader to [T81420] . For an approach in which inter-jet observables are resummedd by nonlinear 
evolution equation which has an analogous structure as in the BK equation we refer the reader 
to [2T]|221I21] and paper where the exact mapping was found [23]. 

The paper is organized as follows. In section 2 we introduce the BK evolution equation 
for the dipole amplitude in the momentum space and perform resummation of unresolved real 
emissions and virtual emissions arriving at a new representation of the BK equation. In section 
3 we present the main result of this paper which is a new evolution equation using the resummed 
BK equation i.e. - the CCFM equation extended by a nonlinear term. 



2 Exclusive form of the Balitsky-Kovchegov equation 



At the leading order in In 1/x the Balitsky-Kovchegov equation for the dipole amplitude in the 
momentum space is [Sj: 
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the linear term can be linked to the process of creation of gluons while the nonlinar term can 
be linked to fusion of gluons and therefore introduces gluon saturation effects. In order to find 
an exclusive form of the BK equation and define a link to the CCFM equation, in the first step 




Figure 1: Plot explaining meaning of variables in BK and CCFM evolution equations. 



we rewrite it as an integral equation following the KMS framework [5] 
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where the lengths of transverse vectors lying in transversal plane to the collision axis are 
k = |k|, I = |1| (k is a vector sum of transversal momenta of emitted gluons during evolution), 
z = x/x'(see Fig. p]), a s = N c a s / V. The impact parameter dependence b of the dipole 
amplitude in momentum space is assumed to be trivial i.e. in a form of a theta function, 
6(R — b) with R defining the target radius. Therefore we suppress it but it is understood 
implicitly. The unintegrated gluon density obeying the high energy factorization theorem [1] is 
obtained from [25H27] : 
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where the angle independent Laplace operator is given by V 
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In order to arrive at an exclusive form of the BK equation first we reintroduce the angular 
dependence to the linear part of eq. ([2]). We obtain 
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where q = 1 — k. Introducing the resolution scale \x and decomposing the linear part of eq. (jlj) 



into resolved real emission part with q 2 > fi 2 and the unresolved part with q 2 < fi 2 , we obtain 
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A convenient way of performing resummation or exponentiation of virtual and unresolved real 
emissions is provided by the Mellin transform defined as 
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Performing the Mellin transform and using in the unresolved part |k + q| 2 k 2 since q 2 < [i 2 
we obtain 

$(oo,k 2 ) = ^(oo,k 2 ) (8) 

+ - / ^[^,|k + q| 2 )% 2 - / , 2 )] + ^ / ^(oo,k 2 )[9^ 2 -q 2 )-9(k 2 -q 2 )] 
oo J q uo J q 

— — pi 



^ / dyy^ 2 (y,k 2 ) 
w Jo 



where in the nonlinear term we changed the variables x/z — >■ y and we integrated over x what 
gives 1/uo in front of the nonlinear part. After combining the unresolved real and virtual parts 
we obtain 

$(oo,k 2 ) = $°(oo,k 2 ) (9) 

+ ^ / |k + q| 2 )% 2 - /, 2 ) - k 2 ) In % 

oo J nq z oo n z 

f dyy^ 2 {y,k 2 )- 
w Jo 

This can be simplified to 

$(oo,k 2 ) = $°(uj,k 2 ) (10) 
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The inverse transform (j7j) can be computed as follows 
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Calculating the residue, changing variables y x/z and using Ar(z, k, fx) = exp (^—a s In ^ In j^j . 
which is called Regge form factor, we obtain: 
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where we introduced 



Finally we obtain 
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Eq. ( fl6|) is a new form of the BK equation in which the resummed terms in a form of Regge form 
factor are the same for the linear and nonlinear part. This form will serve as a guiding equation 
to generalize the CCFM equation to include nonlinear effects which allow for recombination of 
partons. It will also be useful as a starting point in an attempt to solve nonlinear extension of 
the CCFM equation using iterative and Monte Carlo methods [29] . 



3 Towards nonlinear extension of the CCFM equation 

Motivated by the suggestive form of the BK equation where the Regge form factor appears in 
the nonlinear part, in this section we propose an extension of the CCFM equation to account 
for the nonlinearity, what could be viewed as an extension of the BK equation to the large x 
domain. This extension is of course to some extend arbitrary since at present we do not have 
diagrammatic picture of the triple pomeron vertex J301EI] which would account for coherence. 



Our guiding principle is eq. ffTB]) in which the linear part resembles the CCFM equation. The 
expectation is that at the low x limit the solution of nonlinear extension of the CCFM equation 
should approach the solution of the BK equation. Also they should give similar predictions for 
inclusive observables. 



3.1 CCFM evolution equation 

The CCFM equation sums up gluonic emissions with the condition of strong ordering in emission 
angle which allow for smooth interpolation between BFKL limit when z — > and the DGLAP 
limit (in gluonic channel) when z — > 1. The BFKL or the BK are generally applied for the 
calculation of elastic scattering and total cross-sections while due to the dependence on the 
hard scale (which is expressed as the maximal angle) the CCFM equation allows for studies of 
exclusive observables. It is the following 

A(x,k 2 ,p) =A (x,k 2 ,p) (17) 
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The gluon density obtained from the CCFM equation, usually denoted A(x, k 2 ,p), on a level 
of linear equation has interpretation of the gluon density describing parton with longitudinal 
momentum fraction x, and transverse momentum (squared) k 2 which is probed by a hard 
system at scale p. 

The momentum vector associated with z-th emitted gluon is 

q i = oc i p P J r PiPe + Qti (18) 

from which the rapidity and angle of emitted gluon with respect to incoming parent proton 
(beam direction) can be obtained as 

Vi = ~ Mfc) = \ ^ (-) = In ( J^L) , tan § = M- . (19) 

The variable p in (JT7J) is defined via £ = p 2 /(x 2 s) where |ln(£) is a maximal rapidity which is 
determined by the kinematics of hard scattering, y/s is the total energy of the collision. For 
example, for the electron proton scattering s = (pp + p e ) 2 and k' = \k + (1 — z)q\. Using the 
variables £ the angular ordering can be conveniently expressed as: £ > & > > ... > £i > £o? 
where £o = f 1 with being infrared cut off. The momentum q is the transverse rescaled 
momentum of the real gluon, and is related to q by q = q/(l — z) and q = |q|. 
The Sudakov form factor which screens the 1 — z singularity is given by 



The form factor A ns screens the 1/z singularity, in a similar form as the Regge form factor but 
also accounts for angular ordering: 
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where for the lowest value of zq 2 we use a cut off /i. In the literature |32j different forms of the 
non-Sudakov form factor were considered, where the differences came from the approximations 
taken, the one we are using follows results obtained in [11]. For studies of effects of various 
forms of non-Sudakov form factors we refer the reader to [20] • For our purposes it is enough to 
mention that the nonlinearity which suppresses contributions from gluons with low transversal 
momenta makes the dependence on the non-Sudakov's form subleading. This is because the 
non-Sudakov form factor affects mainly the low k part of the gluon distribution. 



3.2 Nonlinear extension of the CCFM equation 

As it has already been stated the motivation to extend the CCFM to account for nonlinearity 
is to be able to study the impact of saturation of partons on exclusive observables. There are 
indications |T8|,[33j that such effects might be significant in for instance production of charged 
particles at HERA or in forward production of di-jets [13] . 

Below we propose the extension of the CCFM equation to include nonlinearity. In order to avoid 
confusion we introduce new notation for the extended by nonlinearity the CCFM equation. We 
also want to recall that the nonlinear extension of the CCFM changes the interpretation of the 
quantity for which the equation is written. It is not longer high energy factorizable gluon density 
A but should be interpreted as the dipole amplitude in momentum space $, denoted from now 
on by £, which depends additionally on a hard scale p. Below we propose the new equation, 
the peculiar structure of the nonlinear term is motivated by the following requirements: 

• the second argument of the 8 should be k 2 as motivated by the analogy to BK 

• the third argument should reflect locally the angular ordering 
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Similarly as in case of the BK equation in order to obtain high energy factorizable unintegrated 
gluon density one applies: 

A non— linear 
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The nonlinear term in (1221) . apart from allowing for recombination of gluons might be under- 
stood as a way to introduce the decoherence into the emission pattern of gluons. This is because 




the gluon density is build up due to coherent gluon emission and since the nonlinear term comes 
with the negative sign it slows down the growth of gluon density and therefore it introduces 
the decoherence. 

We expect the nonlinear term to be of main importance at low x similarly as in case of the BK 
equation. In this limit it will be of special interest to check whether in this formulation of the 
nonlinear extension of the CCFM equation one obtains an effect of saturation of the saturation 
scale as observed in [20]. In the mentioned paper it has beens stated that extension of the 
CCFM equation which mimics saturation by special boundary conditions generates saturation 
scale which saturates itself due to constraints on the phase space coming from maximal allowed 
scale by the kinematics. This effect is of great importance since it has a consequences for ex- 
ample for imposing a bound on amount of production of entropy from saturated part of gluon 
density as observed in [34]. Another interesting limit is the large x region of the phase space 
where the nonlinear effects are believed to be subleading and can be safely neglected. However, 
since the CCFM extends beyond BFKL asymptotics, with the presented formulation one can 
study onset of saturation in a function of x. 

4 Conclusions and outlook 

In this paper we have studied the high energy factorizable evolution equations. We obtained a 
new form of the BK equation where the unresolved real and virtual contributions in the linear 
part of the kernel are resummed. On top of this, both the linear part and nonlinear pieces 
became folded with the Regge form factor. The obtained representation of the BK equation 
might likely be useful in future Monte Carlo solution. However, in order for it to be useful 
for phenomenological studies of exclusive final states one has to take also into account angular 
ordering of gluon emissions. This we do by extending the CCFM equation by nonlinear term 
following the suggestive form of the BK equation in the exclusive representation. In a future, 
the obtained new equation will be investigated with a focus on possibility to solve it using 
Monte Carlo methods. 
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